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Nonsinusoidal current— phase relation in SFS Josephson junctions 
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Various types of the current-phase relation I{^p) in superconductor-ferromagnet-superconductor (SFS) 
point contacts and planar double-barrier junctions are studied within the quasiclassical theory in the limit of 
thin diffusive ferromagnetic interlayers. The physical mechanisms leading to highly nontrivial /(v^) dependence 
are identified by studying the spectral supercurrent density. These mechanisms are also responsible for the 
O-TT transition in SFS Josephson junctions. 
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The relation between the supercurrent / across a 
Josephson junction and the difference ip between the 
phases of the order parameters in the superconducting 
banks is an important characteristic of the structure. 
The form of I{(p) dependence is essentially used for an- 
alyzing the dynamics of systems containing Josephson 
junctions Studying I{ip) also provides information 
on pairing symmetry in superconductors 1^. 

In structures with tunnel-type conductivity of a 
weak link (SIS) the current-phase relation is sinusoidal, 
I{ip) — IcSimp with Ic > 0, in the whole temperature 
range below the critical temperature. At the same time, 
in point contacts (ScS) and junctions with metallic type 
of conductivity (SNS) strong deviations from the sinu- 
soidal form take place at low temperatures T [|| with 
the maximum of achieved at 7r/2 < ipmax < 

The situation drastically changes if there is magne- 
toactive material in the region of weak link. The transi- 
tion from the 0-statc [I{(p) — Ic sin ip with Ic > 0] to the 
TT-statc (Ic < 0) in junctions containing ferromagnets 
has been theoretically predicted in a variety of Joseph- 
son structures ^, |[ ^, 0, ||, ||, |lC , 11 , 12j and experimen- 
tally observed in SFS and SIFS junctions (l|, |l|]. In 
the general case modifications of J((^) do not reduce to 
the O-TT transition. It was shown that tunneling across 
ferromagnetic insulator (Fi) in clean SFjS junctions ||l^ 
or across a magnetically active interface between two su- 
perconductors ll^ may result in a nonsinusoidal shape 
oiI{(p) due to shift of Andreev bound states. Similar sit- 
uation occurs in long SFS junctions with ideally trans- 
parent interfaces in the clean and diffusive [Q, ^ 
regimes. However, in the latter case the effects take 
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place only in a narrow interval of very low tempera- 
tures (due to smallness of the Thouless energy), while 
here we shall consider short-length structures where the 
effects are more pronounced and exist practically in the 
whole temperature range (the role of temperature will 
be discussed elsewhere |]l8|). 

In this letter we investigate anomalies of the I{lp) re- 
lation in several types of SFS structures which allow an- 
alytical solution while have not been fully explored yet: 
the SFcFS point contact with clean or diffusive constric- 
tion as a weak link, and the double-barrier SIFIS junc- 
tion; the ferromagnetic layers are assumed to be thin, 
and the magnetization is homogeneous throughout the 
F part of the system. In particular, we show that the 
maximum of I{ip) can shift from tt/2 < (ySniax < tt to 
< (/3max < 7r/2 as a function of the exchange field 
in the ferromagnet. Previously, current -phase relation 
of this type was theoretically predicted either if super- 
conductivity in the S-electrodes was suppressed by the 
supercurrent in the SNS structure ||l^, |2C , 2l[ or in the 
vicinity of T = in long SFS junctions 17, 81. 

The outline of the paper is as follows. We start with 
studying the SFcFS structure composed of two SF sand- 
wiches linked by a clean Sharvin constriction with arbi- 
trary transparency D. We show that the energy-phase 
relation of this junction can have two minima: at 1^9 = 
and if = TT (while the energy of the junction in the pure 
0- or TT-state has a single minimum — at = or 
</? = TT, respectively). As a result, I{(p) dependence can 
intersect zero not only at ip = and ip = tt but also 
at an arbitrary value ipo from the interval < ipo < tt. 
The salient effects which occur in junctions with clean 
constriction survive averaging over the distribution of 
transmission eigenvalues and thus occur also in diffusive 
point contacts. Physically, the properties of SFS struc- 
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tures are explained by splitting of Andreev levels due 
to the exchange field; to demonstrate this, we study the 
spectral supercurrent. Finally, we show that the same 
mechanism provides shifting of the I{'p) maximum to 
if < 7r/2 in the double-barrier SIFIS junctions which 
can be more easily realized in experiment. 

SFcFS with clean constriction. We start with 
a model structure composed of two superconducting SF 
bilayers connected by a clean constriction with trans- 
parency D (the size of the constriction a is much smaller 
than the mean free path a <C 0- We assume that the 
S-layers are bulk and that the dirty limit conditions are 
fulfilled in the S- and F-metals. For simplicity we also 
assume that the parameters of the SF interfaces 7 and 
7b obey the condition 

7 < max(l,7i3), (1) 

7-B = RbAb/pf^f, 7 = PS^s/pF^F, 

where Rb and Ab are the resistance and the area of the 
SF interfaces; ps{F) is the resistivity of the S (F) mate- 
rial, and the coherence lengths are related to the diffu- 
sion constants Ds{f) as ^s{f) ~ ■\/I?5(F)/27rTc, where 
Tc is the critical temperature of the S-material. We shall 
consider symmetric structure and restrict ourselves to 
the limit when the thickness of the F-layers is small: 



dp < min [ip, \JDfI'2.H^ , 



(2) 



where H is the exchange energy in the F-layers. 

Under condition (|l|) , we can neglect the suppression 
of superconductivity in the S-electrodes by the supercur- 
rent and the proximity effect, and reduce the problem 
to solving the Usadel equations ||2^ in the F-layers 



d_ 
dx 



Gi 



d_ 
dx 



-Gf^f 



0, 



(3) 



with the boundary conditions at the SF interfaces {x 
Tdp) in the form |2|] 

UJ LU 



, ^fGf d ^ ^ 
UJ ax 



(4) 



Gs ^Lo/^Jiu'^ + Al $5(Tdf) = Aoexp(Tiv'/2). 

In the above equations the x axis is perpendicular 
to the interfaces with the origin at the constriction; 
UJ — 7TT{2n+l) are Matsubara frequencies; uj — uj + iH; 
and Aq is the absolute value of the pair potential in the 
superconductors. The function $ parameterizes the Us- 
adel functions G, F, and F: 



Gf{uj) 



(5) 



FfIuj) 



<I>ir(a;) 



FfIuj) = F^{-uj). 



Under condition (g), the spatial gradients in the F- 
layers arising due to the proximity effect and current are 
small. Then we can expand the solution of Eqs. (||)-(^) 
up to the second order in small gradients, arriving at 



$Fi,_F2 = *o exp(=Fi¥'/2), $0 Aot^/VF, (6) 



where 



W = uj + W7sMf^, n=^uj'^ + Al/nT,, (7) 

7-BAf — jBdp/^F, 

and the indices 1 and 2 refer to the left- and right-hand 
side of the constriction, respectively. 

The supercurrent in constriction geometry is given 
by the general expression pj| 



- FiF2)/2 



/^i^ImV 

cRn f^^2^D[l- G1G2 - (F1F2 -f FiF2)/2] ' 

(8) 

where Rf^ is the normal-state resistance of the junction. 
Inserting Eq. (H) in this expression we obtain 



A§ sin Lp 



t^^ -f A2 [1 - sin2(^/2)] ■ 
Finally, the current -phase relation takes the form 
27rr AAl sin Lp 



E 



(9) 



(10) 



A = Al[l-D sin^ ((^/2)] - {-fBM^f 

W{l + lBM^f , 
B = 2ujH-fBM^ (1 + IBM^) ■ 

At small UJ the function A [and thus I{p)] 
changes its sign at finite phase difference 
ipc — 2 arcsin ^[1 — (jBAih)"^]/ D if the exchange 
field is in the range 1 — D < {jBAih)'^ < 1; here h is the 
normalized exchange field, h = H/ttTc- The results for 
I{ip) are shown in Figs.QH and can be understood by 
considering the spectral supercurrent density Im J(£). 
The latter is obtained by the analytical continuation 
in Eq. (||) and is given by a sum of delta-functions 
S{£ — Eb) where Eb are energies of the Andreev 
bound states. At ^bm = the well-known result 



Eb = ±Aoy 1 — -Dsin^(iy9/2) is reproduced, while at 
finite 7_BM the exchange field splits each bound state 
into two (see inset in Figj^). At p — p^ one of these 
split (positive) peaks crosses zero leaving the domain 
e > 0, and simultaneously a negative peak moves from 
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Fig.|^. Current-phase relation in clean SFcFS junc- 
tion with ideally transparent constriction [D — 1) at 
T/Tc = 0.01, 7b M = 1 for different values of the nor- 
malized exchange field h. Inset: spectral supercurrent 
density ai (p = 2n /3 for h — (solid line) and h — 0.4 
(dashed line). 
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Fig.^. Current-phase relation in clean SFcFS junction 
at T/Tc = 0.01, 7flM = 1, /i = 0.8 for different values 
of the barrier transparency D. 



the region e < into e > reversing the sign of the 
supercurrent. 

The sign-reversal of the supercurrent (the O-tt tran- 
sition) can also be achieved at fixed H due to nonequi- 
librium population of levels. This phenomenon has been 
studied in long diffusive SNS (25[ |2|, |2| and SFS junc- 
tions |, |. 

SFcFS with diffusive constriction. To get 

the I{(p) relation for the diffusive point contact [I <C 
a < min(^i?, ^Df/2H)] we integrate p[D)I{D)dD, 
where I{D) is given by Eq. (||) for the clean case 
(note that Rn oc in this equation) and p{D) is 
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Fig.^ Current -phase relation in diffusive SFcFS point 
contact at T/Tc = 0.01, 'Jbm = 1 for different values of 
the exchange field h. 



Dorokhov's density function p{D) 
Finally, we arrive at the result 



AttT 



Re^ 



Aq cos((^/2) 



X arctan 



Ao sin((p/2) \ 



y/W^+Al cos2((p/2) 



(11) 



This expression coincides with the direct solution of the 
Usadel equations, and at jbm = it reproduces the 
Kulik-Omelyanchuk formula for the diffusive ScS con- 
striction p9[ . 

Calculation of /(^s) using the above expression yields 
results similar to those for the clean point contact, how- 
ever the transition from 0- to 7r-state becomes less sharp 
(see Fig.|). 

Temperature dependence of the critical current in 
this case shows thermally-induced O-tt crossover with 
nonzero critical current at the transition point, in agree- 
ment with results of Refs. M, ^ (results for Ic{T) will 
be presented elsewhere This is a natural result 

since the barrier transparency is high and the current- 
phase relation is strongly nonsinusoidal. We note that 



in Ref. 1 13 the measured critical current vanished at 
the O-TT transition point because of the low transparency 
regime (and hence sinusoidal current -phase relation) re- 
alized in that experiment. 

SIFIS. Now we turn to a double-barrier SIFIS junc- 
tion (I denotes an insulating barrier) — this structure is 
easier for experimental implementation than an SFcFS 
junction. In the case of SIFIS, due to dephasing effects 
(this situation is similar to the SINIS junction |3^) the 
supercurrent can not be derived by integrating over the 
corresponding transmission distribution (except for the 
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case of vanishing jbm) and must be calculated by solv- 
ing the Usadel equations. 

We assume that condition is satisfied; then we 
can neglect the suppression of superconductivity in the 
S-electrodes by the supercurrent and the proximity ef- 
fect. In this case the system is described by Eqs. (^- 
(H), although now instead of two F-layers connected 
by a constriction we have a continuous F-layer (at 
—dp < X < dp)- 

We also assume that the F-layer is thin [condition 
(H)] and that 73 dp/^F, hence the spatial gradients 
in the F-layer are small. Then (similarly to the case of 
constriction) we can expand the solution of Eqs. (^)~(^) 
up to the second order in small gradients, arriving at 



<i>F = ^0 cos((p/2) 



Gp — 



. ljGs Aq sin((/3/2) x 
' uiGp 7_B 



6 



(12) 
(13) 



with $0 defined in Eq. (||) [in the final result ( p^ we re- 
tained only the first order in gradients — this accuracy 
is sufficient for calculating the current]. 

Inserting the solution ( |l2| ) , (|l^) into the general ex- 
pression for the supercurrent 



epp 



G 



we obtain 



2ttT 
eR 



'N 



uj>0 



), (14) 



\q sm ip 



A2 cos2(^/2) 
(15) 



(our assumptions imply that Rj\j w 2Rb)- This result 
demonstrates that the SIFIS junction with thin F-layer 
is always in the 0-state.^ Nevertheless, I{ip) is strongly 
modified by finite H (see Fig.^ , especially at small tem- 
peratures. Figure ^ clearly demonstrates that an in- 
crease of H results not only in suppression of the criti- 
cal current, but also in the shift of the maximum 
from (/9max ~ 1.86 at iJ = to the values smaller than 
7r/2. In the hmit of large exchange fields, h ^ 7^1/, 
/((/s) returns to the sinusoidal form. 



^'in the case under discussion when the F-layers are thin and 
the interface parameters obey condition (Q), the phase of the 
pair potential is constant in the S-part and almost constant in 
the F-part, however it jumps at the two SF interfaces ]l2) . The 
two jumps compensate each other in SIFIS with a single F-layer, 
whereas in SFcFS they add up at the weak link thus opening 
possibility for the 7r-state. 
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Fig.^ Current-phase relation in double-barrier SIFIS 
junction at T/Tc — 0.02, "fsAi = 1 for different values 
of the exchange field h. The value ip = 2n/3 will be 
used in Fig.^. 

The physical origin of these results can be clarified in 
the real energy e representation. Making analytical con- 
tinuation in Eq. (^|) by replacement u) — > —i{e + iO), we 
obtain the spectral supercurrent density Im J(e) which 
contains contributions of Andreev bound states with dif- 
ferent energies: 



^ = ^ ^ / i^r) 

Ag sin ip 



Im J(e, H) = Im 



v/Af^v/Agcos2(^/2) 



F, (17) 



e = e + 7BA/(e - H)n(e), ^{e) = J A^ - eV^T,. 



Equation ( |17[ ) implies that at p^ — 2 arccos(7BAf ^) sin- 
gularities in Im J(e) are shifted to the Fermi level. At 
p > LPc the negative singularity in Im J(e) for one spin 
projection crosses the Fermi level and appears in the 
positive energy domain, whereas the positive peak for 
the other projection leaves the domain e > (this pro- 
cess is illustrated in Fig.||). As a result, the contribu- 
tion to the supercurrent from low energies changes its 
sign, and the supercurrent /((/s) becomes suppressed at 
p > Pc (see Fig.^). However, at higher energies e ~ Aq 
modifications in Im J(e) are weak, and the resulting 
I{p) does not change its sign. 

In conclusion, we have studied nonsinusoidal 
current-phase relation in Josephson junctions with thin 
ferromagnetic interlayers and identified the physical 
mechanisms of these effects in terms of splitting of 
Andreev bound states in the junction by the exchange 
field. In particular, we have shown that zero-energy 
crossing of Andreev bound states is responsible for the 
sign-reversal of I{p), which also survives averaging over 
distribution of transmission eigenvalues in the diffusive 
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Fig.^. Spectral supercurrent in diffusive double-barrier 
SIFIS junction with thin ferromagnetic interlayer at 
7flM = 1, f = 27r/3 for two values of the exchange 
field h. The chosen value of <p corresponds to ipc at 
"/BAih — 0.5, and the figure demonstrates that the pos- 
itive peak for one spin projection disappears from while 
the negative peak for the other projection appears in 
the domain e > 0. 

junction. As a result, the energy-phase relation for 
the junction has two minima: dA, ip — and (p — n. 
The phenomena studied in this work may be used 
for engineering cryoelectronic devices manipulating 
spin-polarized electrons and in qubit circuits. 
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